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ABSTRACT 

Context. The quest for hot dust in the central region of debris disks requires high resolution and high dynamic range imaging. Near- 
infrared interferometry is a powerful means to directly detect faint emission from hot grains. 

Aims. We probed the first 3 AU around t Ceti and e Eridani with the CHARA array (Mt Wilson, USA) in order to gauge the 2 fim excess 
flux emanating from possible hot dust grains in the debris disks and to also resolve the stellar photospheres. 

Methods. High precision visibility amplitude measurements were performed with the FLUOR single mode fiber instrument and 
telescope pairs on baselines ranging from 22 to 241 m of projected length. The short baseline observations allow us to disentangle the 
contribution of an extended structure from the photospheric emission, while the long baselines constrain the stellar diameter. 
Results. We have detected a resolved emission around t Cet, corresponding to a spatially integrated, fractional excess flux of 0.98 + 
0.21 X 10"^ with respect to the photospheric flux in the A^'-band. Around 6 Eri, our measurements can exclude a fractional excess of 
greater than 0.6 x 10"- (3cr). We interpret the photometric excess around t Cet as a possible signature of hot grains in the inner debris 
disk and demonstrate that a faint, physical or background, companion can be safely excluded. In addition, we measured both stellar 
angular diameters with an unprecedented accuracy: 0LD(TCet) = 2.015 + 0.011 mas and 0LD(eEri) = 2.126 ± 0.014 mas. 

Key words. Stars: individual: t Cet & e Eri - Stars: fundamental parameters - circumstellar matter - Methods: observational - 
Techniques: interferometric 



1. Introduction 

Planetary systems around main-sequence (MS) stars were first 
indirectly revealed through the detection of their far-infrared 
(IR) photometric excesses by the IRAS and ISO satellites. Early 
interpreted as the signpost of planetary activity, this IR emission 
is thought to arise from short-lived dust grains in gas-free, opti- 
cally thin disks. The production of grains is sustained by aster- 
oid collisions and out-gassing of comets in the first tens of AU 
and by collisions o f Kuiper belt-like bodies at larger distances 
jMever et al.ll20079 . Such grains are also known to form the zo- 
diacal cloud in our own solar system. Given the huge contrast 
with respect to their host star, so far only the brightest nearby de- 
bris disks have been resolved. Single aperture images have high- 
lighted a great diversity of large scale structures comparable to 
our own Kuiper belt (KB) - albeit more massive and extended- 
with a relative void of matter in the central regi ons, where plan- 
ets c ould plausibly have cleared away the dust dSchneider et al.l 
12006 , and recent review by Mann et al. 2006). 

The detections of photometric excesses in the 25- 
100 yum range are much more numerous. Recent surveys, that 
benefit from the high sensitivity of the Spitzer Space Telescope, 
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have revealed not only cold and distant grains in KB analogs, 
but also warmer grains in the inner disk regions where plan- 
ets could be present. For solar-type MS stars, the current pho- 
tometric accuracy of Spitzer places a 3cr detection threshold 
at 100 times the fractional luminosity Ldust/i'* of our zodiacal 
cloud at 70;um, 1000 times at 24 //m and 140 times in the 
8-13 fim range respectivelv. lBrvden et alj (l2006i) have reported 
a disk frequen cy of 13 + 5% at 7 0;ur n, confirming the pr elimi- 
nary studies of lDecin et akl (l2000l) and'Hab ing et all (12001 ). with 
an extended and de eper survey of mature FGK field stars. At 
shorter wavelengths, iBeichman et al.l (l2006l) concluded that as- 
teroid belts 10-30 times more massive than our own are very 
rare, with an 8-13;um excess frequency lower than 2.5% for 
sun-like stars older than 1 Gyr This result confirms the lack of 
warm (Tgr > 300 K) grains tentatively detected at mid-IR wave- 
len gths by previous surv eys. Based on the current sensitivity lim- 
its, |Bryden[etal](|200^ suggest that the statistical distribution of 
positive detections is consistent with most nearby solar-like stars 
harbouring exozodiacal clouds as bright as 0.1 to 10 times our 
own. 

The presence of hot grains in the first AUs of extra-solar 
planetary systems cannot be unequivocally determined by clas- 
sical photometry. The typical accuracy of this method amounts 
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to a few percent at best in rather large fields of view (FOV). 
On the other hand, near-IR interferometers can provide sub-AU 
spatial resolution in a field of view comparable to the size of the 
telescope diffraction pattern. Our observing strategy consists of 
using an optimised set of baseline configurations in order to di- 
rectly measure the excess ratio between the res olved disk and 
the o nly partially resolved stellar photosphere ("D i Folco et al.l 
12004 . This method was recently shown to work wit h the detec- 
tion o f a circumstellar K'-hand excess around Vega dAbsil et al.l 
I2006I) . Fiber-filtered interferometers make it possible to directly 
measure contrasts larger than 100:1 in the near-IR. 

The present paper investigates the case of the brightest sun- 
like stars in our neighborhood, eEridani (K2V, 3.22 pc, K = 1.7) 
and rCeti (G8V, 3.65 pc, K = 1.7). These two stars have very 
close spectral types but different ages: e Eri is younger than 
1 Gyr, while t Get is about twice as old as the Sun (~ lOGyr) 
(ISong et al.ll20"00l:lHabing et al.ll2001blDrFolco et al .''2004). The 
cold regions of their debris disks were imaged at sub-mm wave- 
lengths. The images of e Eri reveal an almost face-on, clumpy 
ring-like structure, peaking at 65 AU (iGreaves et alJI 19981120051 
Poulton et al. 2006). The disk around t Get was detected by 
Greaves et al.l (120041) at 850 fj.m and shows an elongated emis- 



sion, that could be interpreted either as an edge-on disk/ring ex- 
tending out to 55 AU, or as a face-on clumpy structure. The in- 
ferred mas£] amounts to 5-9 Mg around e Eri and 1 .2 around 
T Get, which can be compared with the 0. 1 content of 
cometary-like bodies in our Kuiper belt jGladman et al.ll200rh . 

We propose to further investigate this comparison with the 
inner solar system, by probing the warm dust content of the in- 
terplanetary clouds around these two sun-like stars. This paper 
is the first of a series aiming at directly detecting and charac- 
terising the emission of hot grain s in exozodiacal clouds with 
the GHARA Array interferometer dten Brummelaar eTa n i2005h 
on Mt Wilson (GaUfornia) with the FLUOR bea m combiner 
jGoude du Foresto et al.iri997l: iMerand et an i2006'). The survey 
that we have initiated focuses not only on solar-type stars, as de- 
scribed in the present paper, but also on earlier spectral types, as 
discussed in a forthcoming paper (Absil et al., in preparation). 



2. Interferometric observations and data analysis 

2.1. CHAR A / FLUOR observations 

e Eri and t Get were observed in October 2006 in the TiT'-band 
( 1 .94-2.34 yum). Two pairs of 1 m telescopes were used with sep- 
arations of 34 m (S1-S2 baseline) and 250 m (S2-W1). The short 
baseline allows the large scale emission around the stars (e.g., 
disk, companion) to be resolved, while the long one constrains 
the stellar diameter through the position of the first minimum 
of the visibility curve. The FLUOR beam combining instrument 
measures the squared modulus of the coherence factor (visibil- 
ity) of the two interfering telescope beams, after spatial filter- 
ing of the wavefronts by single-mode fibers. The interferograms 
are recorded through a temporal modulation of the optical path 
difference between the wavefronts with a 250/im scan length 
and a read-out frequency of 500 Hz. The instrumental visibility 
(or transfer function) was estimated using i nterleaved observa- 
tion s of standard stars fr om the catalogues of lBorde et al.l (|2002|) 
and lMerand et al.l ( 12005 ) (see Table[TJ. The interferometric field 
of view is limited to the Airy disk of individual telescopes, i.e. 
about 0.55 arcsec in radius at 2.2 jjm (sUghtly chromatic across 



the ^T'-band). The instrument field of view is further limited by 
the wider FOV of the single-mode fiber, which filters out the 
wavefront corrugations before the beam combination, and can 
thus be affected by a source of incoherent emission. The effec- 
tive FOV results from the overlap integral of the turbulent wave- 
front with the Gaussian-shaped, fundamental mode of the fiber 
(FWHM = 0.8 arcsec). This results in a Hnear FOV of 2.6 AU 
(FWHM) for e Eri and 2.9 AU for t Get. 

2.2. Data analysis and metliodoiogy 

The broad-band squared visibilities C V^) are extracted using 
the FLUOR Data Reduction Software ( Goude du Foresto et al.l 
ll997l:lKervella et al.ll2004t iMerand et al.l l2006) using a wavelet- 
and Fourier-based analysis. The total uncertainty of each visibil- 
ity measurement includes the intrinsic dispersion of individual 

estimations as well as the calibrator induced error on the ab- 
solute visibility. 

The visibilities of the source model are computed as follows: 



,2. „^ / ^fluor(^)[^('^' T,s)/crf^l(B)dcr 



/ 7'fluor(^)[^(^' T.fi)l(Tfd(T 



(T^HA (1) 



' assuming a collisional cascade of km-sized bodies down to 
20 yum grains 



where rFLuoR(o") is the FLUOR chromatic transmission func- 
tion in the /T'-band filter, 'B{cr,T^s)lcr is the Planck function 
(in photons s ' Hz ' m"^) for the stellar effective temperature 
Teff and 'Vj^iB) is the monochromatic visibility of the source at 
the projected baseline B. This formula takes into account the so- 
called bandwidth smearing effect, the result of which is that the 
modeled wide-band visibilities never reach zero. 

While compact photospheres with a typical angular diam- 
eter of Oi, ~ 2 mas are fully resolved with baselines of order 
of 250 m at /I = 2 /im, an extended emission from possible hot 
grains beyond their sublimation distance, over resolved at these 
large baselines, can be resolved with telescope separations as 
short as 10-20m. In order to detect such an IR excess, we first 
compare FLUOR visibilities at the 34 m baseline with the visi- 
bility expected for a purely photospheric emission. To this end, 
we use the high precision angular diameters derived from earlier 
interferometric observations with VLTI/VINGI (Di Folco et al] 
2004) to extrapolate the visibility of the stellar photosphere for 
the 20-30 m range. At such low spatial frequencies, the extrap- 
olated "V^ is insensitive to the photospheric limb-darkening pro- 
file, so the uniform disk approximation can be used. Since the 
stellar surface is almost completely unresolved {'V\ > 0.9), this 
extrapolation can be done with high precision. 

In a second step, a simple model of the source brightness 
distribution is fitted to the full data set, including short and long 
baselines, in order to simultaneously estimate the photospheric 
angular diameter and the contrast of the circumstellar environ- 
ment (GSE), if any. The photospheric limb -darkening follow s 
the, non physical, description proposed by iHestroffen (Il997h : 
h - fj", were fi is the cosine of the azimuth of a surface ele- 
ment of the star where fi - I at the centre of the stellar disk 
and at the limb. This parameterization results in an analytical 
formulation for the monochromatic stellar visibility: 

"V^Ax) = r(v + 1)/^ , with X = nO^Bcr and v = ^ + 1 (2) 

If a significant deficit of visibility is detected at low spatial 
frequencies, the global model takes into account a second com- 
ponent in addition to the star itself. Writing /cse the relative 
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Fig. 1. Squared visibility difference between FLUOR measurements (circles) on the 34 m CHARA baseline and the predicted visi- 
biUty of purely photospheric emission (based on VLTI/VINCI angular diameters). The thin dotted lines correspond to the visibility 
extrapolated from long baseline measurements and the solid lines represent the best fit of a constant "V^ deficit to the new, short 
baseline data (with its 3cr interval - dashed lines). 



brightness of the CSE with respect to the total brightness, the 
final squared visibility can be written as: 



"VHb) = [/cse^cse(B) + (1 - /cse)^*(B)]' 



(3) 



If we assume the CSE to be composed of a uniform emission 
of infinite radius (or at least extended enough to be over resolved 
at the shortest baseline, i.e. VcssiB) - for all B > Bmin), this 
expression can be simplified to: 



T2(B)-(l-/csE)M(fi) 



(B > Bmin), 



(4) 



where 'V^ is the broad-band visibility computed as in ([T]). In 
particular, "V^iB) ^ (1 - 2/cse)^*(B) if /cse <k 1, and the 
visibility deficit A'V^ at short baselines is about twice as large as 
the CSE brightness ratio. 

For the sake of simplicity, and given our limited sampling in 
spatial frequencies, we will assume the source brightness pro- 
files to be circularly symmetric. For e Eri, this is supported by 
the independent estimations of the inclination of its stellar ro- 
tation axis (30 + 3 deg , Cro ll et al. 2006) and of its debris disk 
(about 25 deg, iGreaves et al.ll2005l) . The inclination of r Cet is 
less clear, but both stars have very low rotational velocities, with 
Vrot sin / r espectively equal to 1.4k m.s~' for t Cet and 2.4km.s ' 
for e Eri dValenti & Fischeill2005l) . so that no detectable elonga- 
tion of their stellar surface is expected. 

2.3. Results for t Ceti and e Eridani 

The results of the short baseline measurements are presented 
in Fig. \T\ The projected basehne ranges between 22 and 26 m 
for T Cet and between 24 and 28 m for e Eri. These plots show 
a clear departure from the expected photospheric visibility for 
T Cet. The visibility extrapolated from the earlier estimation of 
the stellar diameter with VLTI/VINCI is plotted as a thin dot- 
ted line for reference (3cr interval). We fit a "V^ deficit (A'V^ - 
"^measured ~ '^predicted-'' which is assumcd to be Constant over the 
projected baseline range: A^^ q^^^^ (-L83 ± 0.35) x 10"^ 
(Xr =0.83). The data do not show any trend associated with 
the baseline azimuth, which supports the assumption of circular 
symmetry. In the case of e Eri, the measurements are consistent 
with the visibility expected from the VLTI/VINCI estimation of 
the stellar angular diameter. Although the dispersion of the mea- 
surements is slightly larger than in the case of r Cet, we can put 



a 3cr upper limit to the visibility deficit in the 20-30 m range 
|A'y2(gEj.i)| < 1.17 X 10-2 ^ 2.4). The rather high value 
is due to the dispersion of the measurements, with a large contri- 
bution from the five points with larger uncertainties but beyond 
2.5cr from the central fit value. The larger dispersion compared 
to T Cet seems to be linked to a degraded coupling of light into 
the instrument fibers, possibly due to poorer atmospheric condi- 
tions. 

The combined data collected at short and long baselines (up 
to 224 m for r Cet and 241 m for e Eri) are then fitted with a limb- 
darkening model of photospheric emission, along with a circum- 
stellar component in the case of t Cet, as explained in details 
in Sect. 12.21 The a values for the limb-darkening (LD) descrip- 
tion of Hestroffe r are de rived from the theoretical 4-parameter 
law tabulated bv IClaret I (12000) for the adopted stellar T^a- and 
metallicity. We find a = 0.15 for r Cet and a = 0.16 for e Eri. 
Because of our sparse sampling of the (m, v) plane, our data do 
not allow us to spatially characterise the shape of the detected 
CSE around t Cet. We therefore consider the simple case of a 
uniformly bright component in the FOV (i.e. a fully resolved 
disk), accounting for a fractional - spatially integrated - flux ra- 
tio /cse- The model visibility, including the bandwidth smearing 
effect, then follows from Eq.lH 



Table 1. Relevant parameters for the calibrator stars. 



Identifier 


Sp. type 




©LD ilo" [mas] 


Target 


HD787 


K5III 


1.80 


2.52 ± 0.03 " 


rCet 


HD 1522 


Kim 


0.92 


3.36 ± 0.04 ^ 


TCet 


HD 6805 


K2III 


0.88 


3.44 ± 0.04 


rCet 


HD 8815 


KOIIIb 


1.48 


2.76 ± 0.03 " 


T Cet , e Eri 


HD 12596 


K2III 


3.12 


1.233 ±0.016'' 


rCet 


HD 12685 


MO/ 1 III 


2.83 


1.450 ±0.019'' 


rCet 


HD 16212 


MOIII 


1.28 


3.11 ±0.03'' 


T Cet , e Eri 


HD 16526 


K4/5III 


2.92 


1.389 ±0.017'' 


rCet 


HD 18071 


G8/K0III 


3.09 


1.056 ±0.014" 


T Cet, e Eri 


HD 29063 


K2III 


2.80 


1.38 ±0.017" 


e Eri 


HD 39853 


K5III 


1.84 


2.38 ± 0.027 ^ 


6 Eri 


HD 220935 


K5III 


2.6 


1.963 ±0.098 " 


rCet 


HD 221745 


K4III 


2.88 


1.426 ±0.019" 


rCet 



" From 
" From 



Borde et all j2002l) 
Merand et all j2005h . 
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Fig. 2. Observations of eEridani (left) and rCeti (right) with FLUOR on CHARA basehnes (S1-S2 and S2-W1). The soUd line is 
the best fit of a source model consisting of a limb-darkened photosphere plus a fully resolved disk. The (Icr) visibility expected 
from the earlier VINCI measurements is represented by the dot-dashed line. The errors on the LD diameters are the Icr statistical 
uncertainties from the fit. The lower part of the panels presents the distance between each single measurement and the best fit 
(dashes stand for the Icr limit and dots for the 3cr) 



For e Eri, the data are consistent with no excess emission 
from the CSE, within a 0.2 x 10"^ uncertainty. The LD model fit 
results in an angular diameter of ©ld - 2. 126+0.009 mas (statis- 
tical error only). The fractional excess emission of its inner de- 
bris disk can be constrained to values lower than 0.6 x 10"^ (3cr 
upper limit). We emphasize that the presence of a faint excess 
due to the CSE not only produces a visibility deficit at short base- 
lines, but also changes the general slope of the visibilities at long 
projected baselines. Therefore it can bias any angular diameter 
estimation by an amount that exceeds the intern al dispersion of 
the m easurement (for more details, please refer to lDi Folco et alj 
l200l . In the case of r Cet, the relative CSE excess amounts to 
/csE = 0.98 + 0.21 X IQ-^ and the de-biased CHARA LD di- 
ameter is estimated to be ©ld - 2.015 + 0.004 mas. Correcting 
the 2004 VLTI data for this bias leads us to the revised value 
©LD = 2.059 ± 0.033 mas. In the next section, we analyse the 
error budget and further compare the CHARA and VLTI mea- 
surements. 

2.4. Discussion 

e Eri and t Cet are very similar in spectral type, magnitude 
and declination, and they also share several calibrator stars, so 
the potential instrumental or calibration bias in our measure- 
ments should be the same for both targets. The detection of a 
short baseline deficit for t Cet, associated with a negative de- 
tection for e Eri, and a comparable statistical uncertainty, rein- 
forces our confidence in the observing method and in the re- 
ality of the detection itself. The observed CSE excess is of 
the order of the excess dete cted aro und the AO dwarf Vega 
(/csE = 1-29 + 0.19 X lO' MAbsil e t al. 2006) and allows the 
stellar angular diameter to be debiased. The derived uncertain- 
ties confirm that we can reach a precision on the CSE excess 
of about 0.2 X 10"^ (Icr) in 10-15 measurements, under median 
atmospheric conditions for K - 2 stars. 

The formal statistical precision of the fitted angular diame- 
ters is very small (5-lOjuas), since we were able to approach the 
first minimum of the visibility function. The total uncertainty 



should also include contributions from the uncertainty of the 
calibrator diameters, the uncertainty of the chromatic transmis- 
sion function Tfluor, and from our a priori assumptions on the 
limb-darkening profiles. For the calibration uncertainty we use 
a conservative value of 0.010 mas. The uncertainty of the stellar 
metallicity (0.2 dex) and Teff (lOOK) leads to a typical uncer- 
tainty for the a-determination of 5 x 10"^, which translates into 
a maximum 2.5 juas error o n the fitted an gular diameters. The 
LD parameters tabulated by Claret (2000) are based on ATLAS 
ID-models of stellar atmospheres. Bigot et al. (2006) have com- 
pared the determination of diameter using the ID ATLAS model 
to those produced by the 3D radiative hydrodynamical simula- 
tions for the star (xCenB, a star whose fundamental parameters 
are very close to those of our targets (KIV, Teff =5260 + 50 K, 
log^ = 4.5, [Fe/H]=0.2). They find that the 3D simulations 
produce a systematically smaller diameter (17 juas difference, or 
0.3 %), mainly impacting the second lobe visibilities. Finally, for 
the uncert ainty of the instrum ental transmission, we use the error 
quoted bv 'Bi got et al. 1 12006) for the VINCI instrument, a con- 
ceptual copy of the FLUOR instrument developed for the VLTI. 
This 0.15 % uncertainty leads to an additional 3 /las error on our 
angular diameters. Quadratically adding the statistical precision 
from the fits and the systematics related to the calibration, to the 
(T-prescription and to the chromatic transmission, we derive a 
total uncertainty of 0.01 1 mas for t Cet and 0.014 mas for e Eri. 
We emphasise that the use of the ID ATLAS models for the a 
priori assumption on the LD profile, which are not constrained 
at all by our observations, may overestimate these diameters by 
an amount comparable to the final uncertainty, compared to 3D 
hydrodynamical models. Measurements in the second lobe of the 
visibility function will be needed to properly estimate the impact 
of the LD profile on the diameter determination. 

A comparison between CHARA/FLUOR and VLTI/VINCI 
measurements, taking into account the total error budget, is pre- 
sented in Table |2] After a correction for the CSE bias, VINCI 
diameter estimations appear to agree with the new values within 
1.3crvinci for T Cet and within O.Scrvinci for e Eri. The preci- 
sion of the Hipparcos parallax contributes 12% to the final 
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Table 2. Global model fits and comparison with earlier estima- 
tions. 



Target 


Instrument 


©LD ±cr [mas] 


/csE±o-(10-^) 




rCet 


CHARA 


2.015+0.011 


0.98 ±0.21 


1.1 




VLTP 


2.078 ±0.031 


<5 






VLTI'^ 


2.059 + 0.033 


0.98 




6 Eri 


CHARA 


2.126 ±0.014 


< 0.6 


2.1 




VLTP' 


2.148 ±0.029 


< 2 





" From lDi Folco et '^( I2004h . 
New fit with fixed /cse value. 



uncertainty on the linear radius of r Cet and 8 % for e Eri: 
RirCet) = 0.790 + 0.005 Rq and R(eEn) = 0.735 + 0.005 Rq. 

3. Nature of the infrared excess 

3.1. An unlikely faint companion around r Cet 

Several physical models can explain the visibility deficit ob- 
served at short baselines. One trivial solution can be an IR point 
source within the FOV of the instrument. Deriving the contri- 
bution of this source from the amplitude of the visibility deficit 
allows us to address the question of its angular separation with 
the scientific target. A maximum upper limit on the source con- 
trast can be set to / ~ A ^^/4 = (0.46 + 0.09) x 10"^, equivalent 
to a maximum observed magnitude niK^ - 7.6 ± 0.3. 

For a physically bound faint source, we can estimate the 
spectral type from the known distance of t Cet (3.65 + . 01 pc) . 
Applying the empirical relations from iDelfosse et alj (l2000l) 
with borderline conditions, it follows that a companion would 
have an absolute magnitude of Mk = 9.7 + 0.3, and thus a 
mass M ^ 0.09 ± 0.01 Mq. Using an age of about lOGyr and 
a metalHcity [M/H] - -0.5. iBaraffe et alj ( [199 8) models indi- 
cate a temperature T^ff = 2800-3000 K and a spectral type M6- 
9V. Orbiting at less than 4 AU, such a massive body would also 
present a clear astrometric and/or radial velocity (RV) signature, 
depending on the inclination of the system. With a mass ratio of 
10 and a maximum separation of 4 AU (period< 7.6 yr), the as- 
trometric signature of the binary would amount to a > 1 10 mas 
(minimum angular semi-major axis). No significant motion was 
detected by Hipparcos a t a 0.8 mas level during the 4-year 
mission jPerrvmanl 1 1 9971) . Moreover, long-term radial velocity 
follow-up of T Cet has reported very stab le measurements over 
about a decade. IWittenmver et alj (|2006|) exclude the presence 
of a Jupiter-mass planet on a 5 AU, circular orbit based on an 
upper limit of the semi-amplitude velocity Kgy ~10-15m.s"'. 
This is also independently confirmed by the CORALIE survey 
which shows a constant radial velocity with a 5m.s"' rms over 
more than 5 yr (Udry, private communication). RV measure- 
ments alone would constrain the orbital inclination (with respect 
to the sky plane) of any 0.09 Mq companion to values smaller 
than 0.21 deg, which is very unUkely and also incompatible with 
the astrometric constraint. 

Furthermore, an IR background source can also be ruled 
out given its low statistical likelihood. The 2MASS survey 
jCutri et al.ll2003h has detected only 355 sources brighter than 
K = 8 in a 5 deg radius patch around t cet. The local surface 
density of such IR sources is thus as low as 3 x 10"^ arcsec"^, 
hence the probability of finding such a faint source in a maxi- 
mum FOV of about 2arcsec^ is 6 x 10"^. The same reasoning 
leads to a probability of 9 x 10"^ for e Eri. In conclusion, al- 



though our interferometric observations alone could in principle 
be reproduced with an additional faint point-like source in the 
FOV, we can confidently exclude the presence of any such com- 
panion. 

3.2. Circumstellar hot dust grains around t Cet 

We show in this section that the detected A''-band excess around 
T Cet can be explained by the emission within our FOV of hot 
dust grains in an exozodiacal disk. Since the spatial frequencies 
probed by the current data set do not allow us to determine the 
exact geometry of the emission, we assume it to be represented 
by a disk of outer radius as large as our FOV, or at least large 
enough to be fully resolved at B ~ 20 m. Following Eq. H] the 
derived integrated flux excess of this inner disk, relative to the 
stellar flux in the /T'-band, is /cse = 0.98 + 0.21 x 10 I Other 
models such as a thin ring or a spherical envelope could also be 
considered and would provide slightly different values for /cse, 
yet have a small impact on our final result. However, the spatial 
characterisation of this CSE is beyond the scope of this paper, 
since it would require many more visibility measurements with 
a continuous sampling of the (m, v) plane. Based on this pho- 
tometric excess val ue, we will use a m odeling approach similar 
to that described in lAbsil et al.l (|2006|) . and successfully used to 
reproduce the resolved emission around Vega. 

We compiled published flux measurements from the visi- 
ble to the sub-mm, resulting in the spectral energy distribution 
(SED) shown in Fig. [3] The s tellar atmosphere SEP is mod- 
eled with a NextGen spectrum (Haus childt et al.iri999 ) with T^fi 
- 5400 K and \ogg - 4.5, and scaled to match the observed V- 
band magnitude (dashed-dotted line on Fig.|3j. The well-known 
far-IR excess beyond 60 /im is also shown in Fig. |3] where it 
has been fitted, for the sake of compariso n, with a modified 
60 K-blackbody following the lGreaves et al.l (120 04) prescription 
(dotted line). This long-wavelength emission component is as- 
sociated with a reservoir of cold material comparable i n size 
to the solar Kuiper belt, and imaged by .Greaves et al.l (|2004|) 
with SCUBA. These cold grains orbit far beyond the 3 AU re- 
gion probed with CHARA, and we will thus concentrate only on 
the warm dust content close to the star. 

The 1-14/im SED of the assumed warm dust population 
(plotted below the star SED on Fig.[3j is obtained by subtracting 
the modeled photospheric emission from the measurements col- 
lected in the literature, except for the jiT'-ban d where we used th e 
excess directly measured by FLUOR. As in lAbsil et all (|2006|) . 
theoretical inner di sk SEDs are ca l culate d using the optically 
thin disk model of lAugereau et al.l (11999) for broad ranges of 
minimum grain radii (amin), dust chemical compositions, disk 
inner rim positions (ro) and surface density power law indexes. 
Compared to the Vega case it appears that the fits are less con- 
strained, in particular because we are lacking accurate 10 yum ob- 
servations, but qualitatively the model requirements are very 
similar We have therefore adopted the grain composition (50% 
glassy olivine, 50% amorphous carbons), the size distribution 
(dn(fl) oc a^^-^da), and the surface density profile (2(r) oc r"'*) 
found to best fit the Vega inner disk. The disk mass within the 
FOV Mdust is then obtained by a least-square fit of the photo- 
metric constraints. The strongest constraints to the fit are the 
A"'-ba nd FLUOR mea surement and the Spitzer/IRS observa- 
tions (Chen et al. 2006, and Spitzer archives data). According 
to Chen et al. (2006), the 10-35//m IRS spectrum is consistent 
with purely photospheric emission within an the uncertainty of 
(5%). The non-detection of a disk at mid-IR wavelengths is here 
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Fig. 3. SED of T Cet. Upper plot: Published photometric measurements are represented by grey squares with Icr uncertainties. 
A NextGen spectrum (dashed-dotted line) is also plotted. The dotted line corresponds to a modified 60 K-blackbody spectrum, 
highlighting the known long-wavelength excesses due to the cold debris disk. Lower plot: residual excess in the l-14//m window 
(diamonds, with Icr error bars), obtained after subtracting the best-fit stellar spectrum. The CHARA ^T'-band direct measurement 
is shown as a black dot. The solid line represents the total emission of the inner disk of hot grains described in Sect. 13.21 while the 
dashed line shows the thermal emission contribution of the disk. In this model, scattered light represents about 20% of the total disk 
emission in the K'-hwd. 



represented as upper limits to the disk emission at four represen- 
tative wavelengths (8, 10, 12 and 14yL/m). 

The best;^'^ values are obtained for minimum grain sizes Omin 
smaller than about 1 /im. Assuming = 0.1 jjm, in order to 
continue the comparison to the Vega case, we find ro - 0.02 AU 
(i.e. 5 stellar radii), the surface density being effectively sharply 
truncated at the sublimation distance rsub(a) for any given grain 
size if ro < rsub(fl)- The fit restricts the possible parameter 
space to ro < 0.1 AU for < I jjm and ro < 0.03 AU 
for flmin > lOyum, although with a lower probability. This sim- 
ple modeling demonstrates that a steep distribution (the density 
scales as r "^) of hot and small grains (flmh, ~ 0.1 //m) can re- 
produce the detected /T'-band excess without generating much 
mid- or far-IR emission. Such a tenuous inner disk would have a 
typical mass of the order of 10"^ (using Omax = 1-5 mm) and 
mean collisional time-scales of a few weeks at the inner disk 
edge (or Mdisk ~ 10"^ Me in the less probable case of grains 
with amin > 10/im). 

For comparison, our zodiacal cloud is estimated from its 
infrared emission to have a total mass between 10"^ and 
10"^ Me dLeinert et al.lll996l: iFixsen & Dwekll2002l) . The sur- 
prising similarity in the orders of magnitude partly comes from 
the much steeper radial dependence of the surfa ce density S(r) oc 
r as compared to r "-^"^ in the Solar system (Kelsall et al.ll 19981 
after integration along the vertical axis). We note however that 
adopting a flatter surface density profile for t Cet (oc r ') in- 
creases the computed mass within the FOV by a factor of only 5. 
The model proposed for t Cet is obviously not a unique solution, 
but it does demonstrate that exozodiacal dust disks around solar- 
type stars might have structures that significantly depart from 
that of the Solar System (see discussion in Sect.|5]l. 



3.3. New constraints on the disk around e Eri 



T Cet, this value translates directly into a maximum fractional 
photometric excess /cse = 0.6 x 10"^ (3cr). Given that the 
stars are very similar in terms of fundamental parameters, we 
have estimated the maximum mass of grains in the inner disk 
in the framework of the model computed for t Cet. According 
to archive Spitz er IRS observations (program ID: 90 see also 
iMarengo Il2005h PI it appears that no significant mid-IR excess 
is detectable below ~ 15;um around e Eri, meaning that the 
photometric constraints for its inner disk emission are similar 
to those of T Cet. Considering that the best-fit model for t Cet 
produces a fractional excess of 0.8 X 10 - at 2yum, a disk mass 
of 10"^ Me (for the steep radial distribution) can be regarded as 
a conservative upper limit for hot grains around e Eri. A tenta- 
tive fit using /cse - 0.35 + 0.25 x 10"^ leads to a preliminary 
value of Mdisk = 5 x 10" '''Me for a r"^ density distribution and 
flmin < 1 yum (10^^ Me for flniin < lOyum respectively). 

We note that the absence of a detectable IR excess, hence 
of hot grains - if confirmed - could be linked to the pres- 
ence of a massive planet (~ 1.5Mjup) orbiting around e Eri 
on an eccentric orb it (a sin/ = 3.4 + 0.4 AU, e - 0.70 + 0.04, 
[Benedict et al.ll2006i) . This giant planet should influence the dy- 
namics of hot grains in the inner region of the debris disk, 
where i t could clear off the dust. Numerical simulations by 
Thebaul t et al. I (I2002h. based on the early planetary parameters 
derived bv iHatzes et al.l (l2000l) . have shown that the region in- 
side the planetary orbit is very hostile to planetesimal accretion 
and terrestrial planet formation. Beyond 0.9 AU, test particles 
are rapidly ejected from the system (this limit might even be de- 
creased d own to 0.6 AU with the larger mass and eccentricity 
revised bv lBenedict et al.ll2p0 6). This is also consistent with the 
recent results oflJones et al.l (12006 ) for the stability of Earth-like 
planets in the habitable zone around e Eri. Inside this dynami- 
cally unstable zone, the relative orbital velocities between plan- 
etesimals are so large (15-40 km/s, Thebault, private communi- 



In Sect. 12. 31 we placed an upper limit on a possible A''-band vis- 
ibility deficit at short baselines of 1.2 x 10"^ for e Eri. Based 
on the same simplistic assumption for the disk geometry as for 



' online presentation available at: 
http://www.stsci.edU/institute/ce nter/information/streaming/archive/l 
NRDD2005/workshopOverview 
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cation) that most collisions would result in disruption rather than 
accretion. Refractory dust grains with collision velocities above 
a thresho ld of 20 km. s~' ma y even be partially or completely va- 
porised (iTielens et al.ll 19941) . The presence of an inner reservoir 
of large bodies to sustain the presence of hot grains in this sys- 
tem seems therefore very unlikely. 

4. Discussion 

4.1. Small grains near the star 

The small hot grains invoked to reproduce the detected emission 
around r Cet are also observed in the inner part of our own plan- 
etary system. In the case of dust located near the Sun, grain prop- 
erties that are relevant to our study can be found in iMann et all 
(12004 . based on in situ spacecraft measurements, remote obser- 
vations at visible and IR wavelength and theoretical modeling. 
Interplanetary dust particles beyond about 0.5 AU are known to 
be concentrated in the ecliptic plane, where they form the zo- 
diacal cloud. Their presence is also inferred down to a few so- 
lar radii by observations of the so-called F-corona during solar 
eclipses. Zodiacal light observations can trace the 1-100 //m size 
particles down to 0.3 AU. They have revealed a flat distribution 
of the dust, with a number density decreasing proportionally to 
about r"'. Inward of 0.1 AU (20 /?o), large changes of the grains 
albedo and polarization in the F-corona suggest a change of dust 
size and composition, possibly linked to the sublimation of var- 
ious dust species. Temporary detections of an IR brightness fea- 
ture during eclipses can also be linked to a significant amount of 
thermal emission, dominated by efficient absorbers like carbon 
grains. These past detections have opened up a debate about the 
possible existence of a dust ring composed of particles of sizes 
a < lOyum, wit h a local density en hancement of a factor 4 in the 
I^Rq region dKimura et al ][T998h . There are therefore similari- 
ties with these models and the grains required to model the inner 
disk around r Cet, however, the detected emission around this 
latter appears much more prominent than in the solar system. 

4.2. Origin of the dust 

Dust supply in the solar zodiacal cloud is supposed to be sus- 
tained by asteroid collisions and evaporation of short-period 
comets, although the respective contributions are unknown. 
Mechanisms invoked to produce the warm dust (Tgr > 300 K) 
observed in the handful of presently detected systems are com- 
prised of: a steady-state evolution of a massive asteroid belt un- 
dergoing collisions and eventually sublimation, the subl imation 
of a " supercomet" captured in a close-in orbit (Beichm an et alj 
l2005h . the recent collision between large planetesimals, and 
the sublimation of a swarm of comets/planetesimals possi- 
bly triggere d by a late heavy bombardment (LHB)-like event. 
I Wvatt et aU (l2007h have investigated these various scenarios 
for the few stars known from Spitzer spectroscopic surveys to 
harbour such warm grains. These are detected through their 
lOyLfm silicate feature, but for the vast majority they are not 
yet not spatially resolved. They conclude that recent and tran- 
sient events that would have caused a population of planetesi- 
mals to be scattered in the inner disk from an outer belt are the 
most likely explanation. Since an LHB is believed to be trig- 
gered when migrating giant planets cross their mutual mean mo- 
tion resonance (1:2 resonance in the case of the solar system, 
[Gomes et al. 2005), this could in principle occur at any time in 
the system's history, depending on the initial conditions for the 
planets relative positions and the planetesimal disk. With an age 



of about lOGyr and a known reservoir of potential comets (one 
or two orders of magnitudes more massive than our KB), t Cet 
is likely undergoing a transient event of dust production, whose 
nature remains to be determined. Precious hints for the origin of 
the dust and its production mechanism might come from future 
detections of long-period giant planets or the detection of silicate 
features. 

4.3. Prospect for DARWIN target selection 

The bright nature of exozodiacal clouds is a potential problem 
for the direct detection of terrestrial planets around nearby stars. 
Future missions like DARWIN or TPF-I/-C, which aim to de- 
tect Earth-like planets around mature solar type stars, will be 
based on nulling interferometry or chronography. The presence 
of small grains in exozodiacal clouds could be a dominant noise 
source in both methods, depending not only on the level of the 
emission but also on its spatial distribution (radial extension, 
disk inclination, presence of structures,...). Earth-like planet de- 
tection in the presence of dust clouds 10 times brighter than our 
own zodiacal cloud might be very difficult, a level still two or- 
ders of magnitude below the current sensitivity limit of Spitzer 
around lOjum. Studying the emission of inner hot dust disks in 
the near- and mid-IR with an increased sensitivity is therefore a 
crucial pre-requisite to the selection of suitable targets for direct 
planet detection. In the case of t Cet, further investigations of 
its inner debris disk are necessary to characterise its spatial ex- 
tension in the terrestrial planets zone and to compare its mid-IR 
brightness to that of the Solar system. Correlating the presence 
of such near-IR emission with excesses at longer wavelengths 
could potentially provide valuable diagnostics for the target se- 
lection process of these future missions. The extension of our 
program to a broader sample of nearby Solar-type stars, either 
with the CHARA/FLUOR facihty or with the VLTI instrumen- 
tation, should make it possible to carry out such analysis in the 
next few years. 



5. Conclusion 

We have used the short and long baselines of the CHARA Array 
and the FLUOR interferometric instrument to probe the close 
environment of two sun-like stars e Eri and t Cet. Our obser- 
vations put a stringent upper limit on the 2 fim emission of pos- 
sible hot dust grains around e Eri and we derive a more pre- 
cise photospheric diameter estimation compared to our previ- 
ous VLTI/VINCI results. Around r Cet, we detect at small short 
baseline a visibility deficit, from which we derive a faint pho- 
tometric excess originating from the circumstellar environment. 
We showed in Sect. 13. II that the case of a bound or background 
companion can be excluded with a high probability, while a sim- 
ple distribution of warm dust grains could be consistent with the 
photometric constraints from the near- to mid-IR. Preliminary 
modeling results show that a grain distribution with micronic or 
sub-micronic minimum sizes, and with surface density peaking 
at a few stellar radii, is best able to reproduce the detected near- 
IR excess together with the lack of significant 10-15/vm emis- 
sion (Spitzer constraint). The 2/vm excess detection calls for 
more efforts to probe the possible lOyum emission of the sus- 
pected grains around r Cet either with MIDI at the VLTI or with 
the nulling mode of the Keck interferometer. Characterising the 
spatial extension of inner exozodiacal disk emissions is a crucial 
precursor science on the road to future terrestrial planet imaging 
missions like DARWIN or TPF 
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